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We report on magnetotransport measurements of nanoconstricted (Ga,Mn)As devices showing 
very large resistance changes that can be controlled by both an electric and a magnetic field. Based 
on the bias voltage and temperature dependent measurements down to the millikelvin range we 
compare the models currently used to describe transport through (Ga,Mn)As nanoconstrictions. 
We provide an explanation for the observed spin- valve like behavior during a magnetic field sweep 
by means of the magnetization configurations in the device. Furthermore, we prove that Coulomb- 
blockade plays a decisive role for the transport mechanism and show that modeling the constriction 
as a granular metal describes the temperature and bias dependence of the conductance correctly 
and allows to estimate the number of participating islands located in the constriction. 
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Narrow constrictions in thin stripes of the ferromag- 
netic semiconductor (Ga,Mn)As display huge magnetore- 
sistance (MR) effects [Il[ai3|4l[5l[^, applicable, e.g., for 
sensors or non- volatile memory elements [5]. The un- 
derlying mechanism causing these MR effects is still not 
unambiguously resolved. Initially it was believed that 
coUinear alignment of the magnetization on both sides of 
the constriction together with a tunnel barrier, formed by 
side wall depletion in the narrow, gives rise to a tunnel- 
ing magnetoresistance effect (TMR) [T]. In these exper- 
iments the in-plane magnetic field was aligned along the 
axis of the stripes. After the finding that the resistance as 
a function of direction and strength of an in-plane mag- 
netic field is reminiscent of the anisotropic magnetoresis- 
tance (AMR), the observed MR effects were ascribed to a 
tunneling anisotropic magnetoresistance effect (TAMR) 
[5], originally observed in an Au/AlOx/ (Ga,Mn)As sand- 
wich structure [7]. There, the angular dependence arises 
from spin-orbit coupling which results in a magnetiza- 
tion dependent density of states (DOS) However, in 
the case of lateral structures, the variation in the DOS is 
^50% [2] and much too small to explain the MR effects, 
orders of magnitude larger [31 U] . A possible explanation 
was based on the assumption that a metal insulator tran- 
sition (MIT) occurs [5j[8], again driven by spin-orbit in- 
teraction and a magnetization dependent overlap of hole 
wavef unctions. On the other hand, experiments on a 
narrow (Ga,Mn)As channel revealed a Coulomb-blockade 
anisotropic magnetoresistance effect (CBAMR), where 
the angular dependence of the resistance is ascribed to 
chemical potential anisotropics [6]. Here we revisit the 
problem of transport across a (Ga,Mn)As nanoconstric- 
tion. In order to shed light on the underlying transport 



mechanism, the bias-, temperature- and gate-dependence 
of the conductance is in the focus of the present letter. 
After a phenomenological explanation for the large MR 
effects and the evidence of Coulomb-blockade (CB) we 
regard the constricted region as a granular metal and re- 
sort to a theoretical model to fit the measured bias- and 
temperature dependence of the conductance. 

We fabricated samples with an individual nanocon- 
striction; an electron micrograph of the central region 
of one of the devices is shown in Fig. [T]^a). The device 
consists of a 20 nm thick Gao.95Mno.05As film grown at 
243° C by low temperature molecular beam epitaxy on 
top of 5 nm Alo.7Gao.3As, 3 nm LT GaAs and a (001)- 
GaAs substrate. After annealing the unpatterned sam- 
ple at 200°C for 8.5 h the Curie temperature, deter- 
mined from the temperature dependence of the resistance 
[H], reached approximately 90 K at a carrier density of 
about 1.8 X 10^''cm~'^. The central area including the 
gates were defined by electron beam lithography where 
cross-linked poly-methyl-methacrylate (PMMA) formed 
the etch mask. Chemically assisted ion beam etching was 
then used to define the (Ga,Mn)As structure. Ti/Au con- 
tacts (not shown in Fig. [ija)) were deposited as source- 
drain and gate electrodes. Each of the devices consists 
of a 3 /^m long and 700 nm wide lead separated by a 
~20 nm wide constriction from a 3 /im long and 100 nm 
wide lead. The structure is aligned along the [100]- 
direction, that is close to an easy axis of (Ga,Mn)As films. 
Apart from that, strain relaxation in the stripes turns 
the easy axis towards the leads [TU]. To determine the 
switching fields for a magnetization reversal in the wide 
and narrow stripes via the AMR effect some of the devices 
have additional potential probes on each stripe. Below 
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FIG. 1; (color online) (a) Electron micrograph of the cen- 
tral part of a device, tilted by 40°. (b) Polar plot of Rsd 
at 1.6 K showing the strong anisotropy of Rad for sample A 
as a function of the magnetization direction. The measure- 
ment was done in a high magnetic field of 300 mT. (c) MR 
of sample A for a = Q° and a = 30° at 1.6 K. The switching 
fields of ±28 mT and ±39 mT (for a = 0°) correspond to the 
magnetization reversal of the broad and narrow lead, respec- 
tively Comparing the resistance values of (b) and (c) allows to 
deduce the magnetization alignment in the constriction. The 
configurations for a HR state together with the magnetization 
angle in the constriction is sketched in (d). 



we present data of two samples A and B. While sample A 
has no side gate, sample B has additional side electrodes 
(as presented in Fig. [l|^a)), separated by ~130 nm from 
the constriction, to tune the electrostatic potential in the 
constriction. Magnetotransport measurements were car- 
ried out in a ^He bath cryostat or a '^He/^'He dilution 
refrigerator, each equipped with a superconducting mag- 
net. The angle a between the applied in-plane magnetic 
field and the current direction could be varied by rotat- 
ing the sample inside the cryostat. The measurements 
were done using dc-technique in a two-point configura- 
tion: a constant voltage Vsd between source and drain 
was applied and the resulting current / was measured 
employing a current amplifier. 

First we demonstrate that the MR effects observed in 
(Ga,Mn)As nanoconstrictions can be phenomenologically 
explained by the magnetization alignment in the constric- 
tion. Figure[T]^b) shows a polar plot of the resistance Rsd 
of sample A measured at 1.6 K by rotating the sample in 
an in-plane magnetic field of 300 mT. This field strength 
is large enough to align all magnetization vectors into the 
external field direction. R^d is strongly anisotropic with 
resistance changes of up to a factor of 20. Neither AMR, 
which in bulk (Ga,Mn)As is on the order of a few percent, 
nor TMR (magnetization in the leads is always parallel) 
or TAMR (change of DOS too small) can explain such 



drastic resistance changes [1] . Figure [TJc) displays the 
MR as a function of the in-plane field during a magnetic 
field sweep along the [100]-structure axes (a = 0°) and 
for the direction causing the largest MR effect (a — 30°). 
The switching fields of the spin- valve like signal perfectly 
agree with the magnetization reversal fields of the 700 nm 
and 100 nm wide stripe that have been detected simul- 
taneously by a four-point measurement using a sample 
with additional voltage probes on each lead (not shown 
here). Note that the switching fields are different in the 
wide and narrow stripe due to a different strain relax- 
ation. The large resistance change in Fig. [ijc) can there- 
fore be ascribed to the consecutive magnetization rever- 
sal of the 700 nm and the 100 nm stripe. The maxi- 
mum resistance observed in the high field experiment of 
Fig. [ijb) is in good agreement with the high resistance 
(HR) measured in the experiment of Fig. [ijc). This sug- 
gests that both effects stem from the same origin and that 
Rsd is linked to the magnetization direction in the con- 
striction. Thus, by comparing the resistances of the polar 
plot of Fig. [TJb)with the spin- valve like signal observed 
during a magnetic field sweep we can deduce the mag- 
netization direction in the constriction and explain the 
MR trace by means of the magnetization alignments in 
the device. The magnetization configurations for the HR 
states marked with open circles in Fig. [ijc) are sketched 
in Fig.[ljd). So, at low external magnetic fields, the rela- 
tive alignment of the magnetization in the wide and nar- 
row stripe involves a distinct magnetization orientation 
in the constriction and thus determines the resistance. 
Within this picture the MR is easily explainable for an 
applied magnetic field along the stripe axis (i.e., a = 0°). 
Here, the HR state occurs due to the 180° magnetization 
reversal in the 700 nm stripe at -28 mT whereas the mag- 
netization in the 100 nm stripe still remains in the [100] 
easy axis [see sketch 1 in Fig. [ijd)] until its cocrcitive 
field is reached at -39 mT. The antiparallel alignment of 
the stripes causes the HR due to the associated orien- 
tation of the magnetization in the constriction and not 
due to the TMR effect. But also more complex features 
of the MR trace that appear for magnetic field angles 
a 7^ 0° can be understood. Consider, e.g., the MR in the 
case of a = 30° , shown in Fig. [ijc) . The resistance in- 
crease within the magnetic field range marked with 2 and 
3 arises from a coherent magnetization rotation in the 
constriction due to an increasing magnetic field strength 
along the 210° (= —30°) direction. This is illustrated in 
the cartoons 2 and 3 of Fig. [ijd). The magnetizations 
in the leads remain unchanged in the considered mag- 
netic field window due to the strong uniaxial anisotropy 
of the stripes. Increasing the magnetic field strength fur- 
ther leads to an abrupt 180° magnetization reversal in 
the 100 nm stripe causing the abrupt resistance change 
due to the resulting magnetization alignment in the con- 
striction towards 180°. 

We now address the microscopic origin of the large MR 
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effects and discuss the results with respect to the avail- 
able models, magnetization induced MIT as proposed in 
[5] and Coulomb-blockade assisted as reported in [6]. At 
temperatures below T ~ 30 K the resistance of our pat- 
terned samples increases continuously with decreasing 
temperature (not shown here), both for the high resis- 
tance (HR) and low resistance (LR) state, thus indicat- 
ing that the huge MR effects, which vanish above ~ 30 K, 
only exist in the insulating regime. A picture involving 
a magnetization driven MIT, for which a different tem- 
perature behavior for the HR and LR state would be 
expected, is therefore at odds with our experiment. In 
the CBAMR picture [6], however, the system is expected 
to be insulating both in the HR and LR state and is 
consistent with our results. The model is based on the 
assumption that disorder potential fluctuations together 
with side wall carrier depletion create small, isolated is- 
lands in the constriction. 

To prove that in our (Ga,Mn)As point contacts 
Coulomb- blockade is at work we investigated samples 
with additional side-gate electrodes. The non-linearity of 
the I-V characteristic of device B, plotted in Fig. |2ja), 
clearly depends on the gate voltage. The corresponding 
conductance G [Fig. |2]^b)], measured in an in-plane field 
of 0.73 T [H] along the [100]-direction at T = 0.55 K, 
displays pronounced oscillations, indicating Coulomb- 
blockade. These conductance oscillations have an ir- 
regular spacing, suggesting that more than one island 
is formed in the constriction. From the average oscil- 
lation period AVg in Fig. ^h) we estimate a capaci- 
tance Cg ^ 0.55 aF between gate electrode and islands. 
For a very rough estimate of the island area we use 
the charging energy of some meV taken from the I-V- 
characteristics. This gives a total capacitance Cs of 
order 10 aF. By approximating the island as a sphere 
we thus expect a radius of order 10 nm consistent with 
the constriction size of our device. The pronounced CB 
diamond-like structure in Fig.[2](c) highlights the role of 
CB in transport across (Ga,Mn)As point contacts. The 
anisotropy of the resistance and its dependence on Vg 
is most clearly seen in Fig. ^d) where polar plots of 
Rsd are shown for three gate voltages. With different 
gate voltage the direction of B and hence of the magne- 
tization where the largest resistance is observed changes. 
Hence, both the gate voltage and the magnetization are 
independent parameters which allow us to tune the con- 
ductance through ferromagnetic (Ga,Mn)As islands. The 
underlying physics has been described by Wunderlich et 
al. |B] within the CBAMR model and is due to (different) 
changes in the chemical potential inside and outside of 
conducting (Ga,Mn)As islands. Due to the strong spin- 
orbit coupling in (Ga,Mn)As the chemical potential in the 
island depends on both the gate voltage and the mag- 
netization direction. These effects might be enhanced 
by the magneto-Coulomb effect [TT]. The different re- 
sistance anisotropy for devices A and B at Vg = (see 



a = 


0° 






T = 


0.55 K 






B = 


730 mT 










U = 

g 

U = 

g 


280 mV 
100 mV 



(b) 



0.4 

\ 0-1 
O 0.0 



H 

AV 



IP 



^0° 



ru/l 




330 



-1.0 -0.5 0.0 
V (V) 



FIG. 2: (color online) (a) I-V characteristics of sample B for 
two different gate voltages showing that the device acts as a 
transistor, (b) Dependence of the conductance on the gate 
voltage for a = 0° and Vsd = 3 mV. The oscillating behavior 
is ascribed to a CB transport. AVg stands for the average 
oscillation period, (c) Color plot of the conductance versus 
gate and source-drain voltage resulting in a CB diamond-like 
structure. The diamond pattern changes for different magne- 
tization directions, here shown for 0°. (d) Polar plot of Rsd 
for different gate voltages. All measurements were carried out 
at B = 730 mT and T = 550 mK. 



Fig. [T]^b) and|2jd)) can be explained by a different elec- 
trostatic environment in the vicinity of the constriction 
due to thermal cycling [T^ . 

As a result of the gate measurements we conclude that 
more than one conducting island is located within the 
nanoconstriction. To elucidate the transport mechanism 
further we studied the temperature and bias voltage de- 
pendence of the conductance G across the constriction. 
Corresponding data are shown in Fig. |3] for both, HR 
and LR state. While in Figs. [3] (a) and (b) the conduc- 
tance is plotted as a function of 1/VT, Figs, [sje) and 
(f) display G as function of T for different source-drain 
voltages Vsd- The conductance shows a strong bias de- 
pendence [see also Figs.|3jc) and|3]^d)]. Such behavior is 
expected once Coulomb-blockade is involved in the trans- 
port process. As our system is reminiscent of a granular 
metal film - metal grains separated by insulating shells 
- which have been investigated intensively in the past 
(see, e.g., [TS] and references therein) - we resorted to 
a corresponding model. We went back to the model of 
Abeles et al. [13] and applied it to our system. Due to 
the surface depletion in the nanoconstriction we assume 
that metallic spots of (Ga,Mn)As are separated from each 
other by thin insulating regions, thus justifying the gran- 
ular metal assumption. As shown below the modeled 
bias and temperature dependence of the conductance fits 
perfectly the experimentally observed data. The granular 
model distinguishes two regimes reflecting the interplay 
of the involved energy scales: the thermal energy fc^T, 
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the charging energy Ec = e'^/2Cs and the energy asso- 
ciated with the average voltage drop across one island 
eAV = eVsd/N. Here, iV indicates the number of is- 
lands whereat the voltage drop occurs. For low electric 
fields where eAV is much smaller than fcsT the model 
is mainly based on charge-carrier generation by thermal 
activation and predicts: G{T) = Go exp{— 2(C/Vfcs?")}- 
Here, C is a constant proportional to the charging energy. 
However, our experiments are clearly not in this low-field 
regime since the data do not follow the InG oc T~^/^ 
dependence over the entire temperature range shown in 
Figs.[3];a) and[3];b), both for the HR and LR state. This 
is consistent with a simple estimate: for our highest tem- 
perature of 0.8 K corresponding to ksT/e ~ 0.07 mV, 
and lowest applied bias voltage, Vsd = 2.1 mV, the low 
field regime requires that the number of islands N is much 
larger than 30. Considering the dimensions of our con- 
striction this is highly unlikely. Thus, we have to examine 
the high-field regime {eAV > ksT) where field induced 
tunneling becomes important and increases the carrier 
density. In this case the conductance is given by [T3] 
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with the electron charge e and the constant Vq depend- 
ing on the charging energy. In the limiting case of zero 
temperature the conductance is given by 



G{Vsd) = Goo exp i-Vo/Vsd) 



(2) 



Hence, plotting InG versus I /Vsd at the lowest temper- 
ature allows to estimate both Goo and Vq directly from 
experiment. The corresponding data are shown for the 
HR and LR state in Fig.jsjc) and (d), respectively. The 
only undetermined parameter in Eq. [T] is the number of 
islands N which we use as the only free fit parameter. 
The experimental data at temperatures below 1 K, plot- 
ted for different bias voltages and for HR and LR states 
in Figs. [Sje) andjsjf), respectively, can be very well fit- 
ted by using N — 2. Eq. [T] describes the temperature 
dependence for all three traces correctly, both for HR 
and LR state, using within the error margins the same 
values for Goo and Vq given in the caption of Fig. [3j The 
high sensitivity of the fits with respect to the number of 
islands is demonstrated in Figs.jsjg) and (h). Assuming 
that the constriction contains two islands and using the 
previously estimated island size of 20 nm diameter sug- 
gests an individual island diameter of order 10 nm. The 
granular metal model, though at the fringe of applica- 
bility given the low N value, describes surprisingly well 
the bias and temperature dependence of the conductance. 
This provides further evidence that Coulomb-blockade 
plays a prominent role in transport through (Ga,Mn)As 
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FIG. 3: (color online) [(a - b)] logG versus T"^/^ plot 
in a HR and LR state realized by applying a B-field of 
730 mT along 120° and 90°, respectively, showing that our 
experiments are not consistent with the low-field regime of 
the model for granular metal films, [(c - d)] Voltage de- 
pendence of G showing a good agreement with the model 
for the high-field regime. From the linear fit we obtain 
\n[Goc{HR)/10-^e'^/h] = 32 ± 2, Vb = (0.15 ± 0.01)V 
for a = 120° (HR) and In [Goo{LR)/10-'^e^ /h] = 18 ± 2, 
Vo = (0.038 ± 0.003) V for a = 90°. [(e - f)] Experimental 
data and numerical fits for the temperature dependence of G 
for different Vsd for HR and LR state, respectively. is the 
fitting parameter and indicates the number of islands which 
are schematically illustrated in the inset. The sensitivity of 
A'^ is demonstrated in (g) and (h). 



nanoconstrictions and is responsible for the huge magne- 
toresistance effects observed. 
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